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ABSTRACT 

VO ■ The currently feasible method of detection of Earth-mass planets is transit 

photometry, with detection probability decreasing with a planet's distance from 
the star. The existence or otherwise of short-period terrestrial planets will tell 
us much about the planet formation process, and such planets are likely to be 
detected first if they exist. Tidal forces are intense for short-period planets, 
and result in decay of the orbit on a timescale which depends on properties of 
the star as long as the orbit is circular. However, if an eccentric companion 
planet exists, orbital eccentricity (ej, % for inner orbit) is induced and the decay 
timescale depends on properties of the short-period planet, reducing by a factor 

^ of order 10 5 e? if it is terrestrial. Here we examine the influence companion planets 

have on the tidal and dynamical evolution of short-period planets with terrestrial 
structure, and show that the relativistic potential of the star is fundamental to 
their survival. 

Subject headings: relativity — planetary systems — celestial mechanics — stars: 
late-type — stars: low-mass, brown dwarfs 



1. Introduction 

In the past few years, the radial velocity method has been used to detect more than 100 
extrasolar planets (Kron et al. 2000) with minimum masses in the range 0.11 to 17 Jupiter 
masses. With the current lower detection limit of 3 m/s, discovering Earth-mass planets 
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with periods longer than a few days, this way is untenable. In principle, Earth-mass planets 
with periods less than 3 days may be marginally detectable with 1 m/s velocity precision 
(Narayan et al. 2004). But a large number of observations is needed to reduce the noise 
associated with the stellar radial velocity "jitter" (Saar et al. 1998). 

Hopes for detecting Earth-mass planets in the near future lie with transit photometry 
and microlensing. So far, four gas giants have been observed crossing the faces of their stars 
(Brown et al. 2001; Konacki et al. 2003, 2004; Bouchy et al. 2004), three of which were 
discovered using microlensing techniques (Konacki et al. 2003, 2004; Bouchy et al. 2004). 
While the transit of a gas giant results in a maximum reduction of the star's brightness of 
around 1-2%, an Earth-mass planet would dim the star's light by only a few parts in 10,000. 
Nonetheless, such high-precision photometry is currently feasible in the MOST satellite and 
will be employed in the Kepler (Koch et al. 2001) mission. Simple geometrical arguments 
give the probability of seeing a planet cross the face of its parent star along our line of 
sight as R*/2a, where R* is the radius of the star and a is the distance of the planet from 
the star during the transit. Thus it is most likely that short-period Earth-like planets will 
be discovered first if they exist, especially since a weak signal can be strengthened by the 
integration of many orbits. This is not only true for planets discovered directly with transit 
photometry, but also those discovered using microlensing techniques since followup transit 
photometry is necessary for their confirmation and study. 

In Section 2 we review the process of tidal damping, in Section 3 we discuss the dynam- 
ical effect of the presence of a companion, and in Section 4 we present a discussion. 



2. Orbit Evolution due to Tidal Damping 

During the formation epoch, tidal interactions between Jupiter-mass protoplanets and 
their nascent disks can cause them to spiral in toward their host stars (Goldreich and 
Tremaine 1980; Lin et al. 1986), a migration scenario (Lin et al. 1996) adopted to account 
for the origin of extrasolar planets with periods of a few days (Mayor and Queloz 1995). 
Other scenarios have been suggested (Rasio and Ford 1996; Nagasawa et al. 2003; Murray 
et al. 1998), all of which involve significant excitation of orbital eccentricity. Once a planet 
is in the vicinity of the star, however, tidal forces will tend to circularize the orbit as well as 
synchronize the planetary spin (Rasio et al. 1996). Earth-mass terrestrial planets may also 
be brought close to their host stars via one or more of these processes and their discovery (or 
otherwise) will shed much light on the planet formation process. For example, along their 
migration paths, gas giants induce residual terrestrial planets to undergo orbital evolution 
through their resonant and secular interaction. The detection of both short-period gas giants 
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and terrestrial planets around any host stars would provide overwhelming support for the 
core-accretion scenario. 

Short-period planets are also capable of raising a substantial tide on their host stars. 
In order to avoid noise introduced by stellar jitters, the present planet-search campaigns 
focus on target stars with quiet chromosphere and slow spin rates. With one exception 
(Butler et al. 1997), all stars observed hosting short-period ("hot") Jupiters have spin periods 
considerably longer than the orbital period of the planet. Mature solar-type stars tend 
to be sub-synchronous because stellar winds carry away spin angular momentum through 
magnetic braking (Soderblom et al. 1993). Within these slowly spinning stars, the induced 
tidal oscillations act to dissipate their acquired tidal energy, the tidal response tends to lag 
the line of centers of the star and planet, and the resulting torque tends to spin up the star 
at the expense of orbital energy and angular momentum. 

The simplest quantitative description for tides is the equilibrium tide model (Goldreich 
and Soter 1966; Hut 1981; Eggleton et al. 1998) which is adequate for analyzing the tidal 
response of terrestrial planets. In gaseous planets and stars, tidal perturbation can excite 
resonant responses which strongly enhance the dissipation rate (Zahn 1970). In radiative 
stars, the response is mostly in the form of g-mode oscillations whereas in fully convective 
envelopes of rapidly spinning gaseous planets and stars, it is in the form of inertial waves 
(Ogilvie and Lin 2004). Since we are primarily interested in the fate of the terrestrial planets 
in this paper, we adopt here the equilibrium tide model. 

To leading order in the orbital eccentricity, the timescale, r , for orbital decay is given 
by (Mardling and Lin 2002) 
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are the contributions to orbital decay (or expansion) from tidal dissipation within the star and 
planet respectively, and the assumption has been made that the planet spins synchronously 
with the orbital motion. Here P Q rb and P spm are the orbital and stellar spin periods respec- 
tively, <2j is the semi-major axis of the orbit, n,i is the mean motion (orbital frequency), M i5 
Ri, and Q\ are respectively the mass, radius and modified Q-value of the planet, the latter 
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containing information about damping efficiency and rigidity (Goldreich and Soter 1966) 
(see Table 1), and M*, i?*, and Q'^ are the corresponding values for the star. If the orbit is 
perfectly circular the planet makes no contribution to orbital shrinkage and the orbital decay 
timescale depends entirely on properties of the star. Otherwise the planet dominates the 
tidal decay process, and this is particularly true if it is Earth-like. In the case of slow stellar 
rotation and small eccentricity, the ratio of contributions to orbital decay from the planet 
and the star is 12.5ef for a Sun- Jupiter pair, 1.3 x 10 5 ef for a Sun-Earth pair, and 2.2 x 10 5 e^ 
for an M dwarf-Earth pair, where we have used Q'- values from Table 1. The discrepancy 
between the Jupiter system and the Earth systems is mostly due to the factor 10 4 difference 
in Q'-values of gas giants and terrestrials. Thus even a small orbital eccentricity allows a 
planet to dominate the tidal decay process if it has terrestrial structure, and can reduce the 
decay timescale considerably. 

The timescale, r e , for eccentricity damping (or excitation) is given by an expression 
similar to Eqn. 1 (Goldreich and Soter 1966; Mardling and Lin 2002) except that it does not 
depend on e$ to leading order in e*. Nonetheless, unless the star spins rapidly, the planet 
dominates the tidal circularization process (Dobbs-Dixon et al. 2004). Table 2 lists orbital 
decay and eccentricity damping timescales, together with other system parameters including 
transit probabilities for various real and hypothetical single planet systems. Included are 
the three of the shortest period planets discovered to date as well as HD209458. The values 
of a,i for the three hypothetical systems (Sun-Jupiter, Sun-Earth and M dwarf-Earth) were 
chosen to indicate how close to the host star we can ever expect to find planets in these 
types of systems. Note especially the steep dependence of r a on the eccentricity when the 
planet is Earth-like. 



3. Eccentricity Excitation due to a Companion 

HD209458 is included in Table 2 as one of only two transit systems discovered so far 
(Brown et al. 2001), the data for which allows an estimate of the planet radius of 1.35i?j 
which we use here to estimate timescales. The planet's relatively large size can be accounted 
for by inflation due to tidal dissipation of energy (Bodenheimer et al. 2001) if e« ~ 0.03 and 
Q\ = 10 6 , with the corresponding orbital decay timescale being comparable to the age of the 
Universe. The observationally estimated eccentricity is small (< 0.05) but not zero (Fischer, 
private communication), even though r e is smaller than the estimated age of HD209458 
(see Table 2). In order to resolve this paradox, the existence of a second planet has been 
conjectured (Bodenheimer et al. 2001) which secularly excites the eccentricity (Murray and 
Dermott 1999). 
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This eccentricity-excitation mechanism puts severe constraints on the existence of sys- 
tems containing short-period terrestrial planets with companion planets. Figure la shows the 
evolution of the eccentricity of an Earth-mass planet in a 2.3 day orbit (aj=0.02 AU) around 
an M dwarf of mass 0.2M Q which spins with a period of 40 days. It has a Jupiter-mass 
companion at a distance 0.7 AU which itself has an eccentricity of 0.2. The initial eccentric- 
ity of the inner planet is zero, however, it is excited by the presence of the companion to a 
maximum value which depends on the strength of additional perturbing accelerations which 
include the spin and tidal bulges of the star and planet, as well as the relativistic potential 
of the star (Mardling and Lin 2002). 

In the absence of any of these accelerations, the inner eccentricity is excited to a value of 
0.015 (the black curve). Taking half this value to represent the equilibrium inner eccentricity 
(see following paragraphs), the orbital decay timescale would be 7.1 Gyr. However, if all the 
perturbing accelerations listed above are included, the maximum inner eccentricity is only 
0.0035, corresponding to an orbital decay timescale of 84.3 Gyr. This suggests that relativity 
plays a major role in the survival of short-period terrestrial planets with companions (given 
that Nature produces such systems in the first place). In fact it is the contribution it makes 
to the apsidal advance of the inner orbit which is responsible for the inhibition of the growth 
of the inner eccentricity (Holman et al. 1997). The fractional contribution the relativistic 
potential makes to apsidal advance compared to that made by a companion planet is given 
by 7/(1 + 7), where 7 = 4Z 2 (M*/M )(a /aj) 3 /(l — ef), with M the mass of the companion 
planet and Z = aifii/c the ratio of the orbital speed of the inner orbit to the speed of light 
(Mardling and Lin 2002). For the case illustrated in Figure la, the relativistic potential is 
responsible for 95% of the apsidal advance of the inner orbit, while in the case of Mercury's 
orbit around the Sun, it contributes only 7%. 

Even more extreme is the case in which a short-period terrestrial planet has an Earth- 
mass companion. For example, for a Sun- Earth- Earth system with semi-major axes 0.03 AU 
and 0.5 AU, the relativistic potential is responsible for 99.9% of the apsidal advance of the 
inner orbit. If the eccentricity of the outer planet is 0.3 the orbital decay timescale is 702 
Gyr, while in the absence of relativity it would be only 0.47 Gyr. 

Using equations which govern the secular evolution of the orbital elements for a dissi- 
pationless point-mass coplanar system (Mardling and Lin 2002; Murray and Dermott 1999; 
Wu and Goldreich 2002) and assuming small values for the eccentricities and ai/a one can 
obtain an estimate for the variation of the inner eccentricity, be^ which also holds when the 
minimum eccentricity is zero at which point e, is discontinuous as in Figure la. This esti- 
mate will vary slightly for moderately non-coplanar systems, and does not apply to resonant 
systems (Murray and Dermott 1999; Novak et al. 2003). If tidal dissipation is taken into 
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account, the system behaves as a damped autonomous system with the familiar circulatory 
and libratory behaviour (Fig. lb). For a given set of initial conditions the system evolves to 
a fixed point corresponding to Wi — w = 2nn, where n is some integer and wi and w are the 
longitudes of periastron of the inner and outer orbits respectively, and a finite equilibrium 
inner eccentricity, e\ , given approximately by 

e ( e( i) = s( a »/ a °) e ° . = 8e-l2 (4) 

|1 -(M,/M ) v^K + 7 | 

Note that as well as the relativistic potential, 7 may include contributions from other per- 
turbing accelerations such as tidal and spin bulges. Figure lc shows the dependence of ef 
on Mi/M in the absence of relativistic and other perturbing accelerations (7 = 0). The solid 
curves were obtained by integrating the Newtonian equations of motion for point masses (av- 
eraged over the inner orbit), while the dashed curves are given by Eqn. 4. Agreement is 
good except near the point corresponding to Mi/M = ^Ja /a,i. The inclusion of relativistic 
effects introduces an additional scale which is proportional to rrii/ai, and this is illustrated 
in Fig. Id for the case a^/ai = 10. Only in one case is the discrepancy between the full 
solution and that given by Eqn. 4 evident. 

Figure 2 consists of grey-scale plots of the orbital decay timescale of an Earth-mass 
planet orbiting a solar-mass star at 0.03 AU, as a function of the eccentricity, e , and semi 
major axis, a , of a Jupiter-mass perturber. The equilibrium eccentricity is calculated using 
the full governing equations (Mardling and Lin 2002), and this is then used in Eqn. 1 to 
obtain the decay timescale. The top panel is for the case where no perturbing accelerations 
are included, while the bottom panel includes the relativistic potential of the star. It is clear 
that relativity allows many systems to survive which would otherwise have suffered tidal 
destruction during the current lifetime of the star. 



4. Discussion 

Of particular interest amongst short-period terrestrial systems will be those with low- 
mass stars for which the habitable zone (planet surface temperature in the range 0-100°) is 
at a distance where the transit probability is not negligible. Table 2 lists several hypothetical 
systems composed of a O.2M M dwarf and an Earth-mass planet for which the habitable 
zone is around 0.04 AU. Planets in such close proximity to their host stars will be tidally 
locked so that one side of the planet is never directly heated. However, systems for which the 
habitable zone is further out may have planets locked in other spin-orbit resonances. Mercury 
is locked in a 3:2 spin-orbit resonance which relies on its permanent slight departure from 
sphericity as well as its substantially non-circular orbit (Goldreich and Peale 1966). Such 
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planets would be heated more evenly and hence be perhaps better candidates for detecting 
the signatures of life in their atmospheres. 

While the Q-values of putative short-period terrestrial-type planets may turn out to be 
somewhat higher than the Earth's, either because of the absence of oceans or because of the 
different temperature profiles caused by the close proximity to the parent star, it is clear that 
a star's general relativistic potential must play a major role in the survival of such planets, 
and hence it is vital that it be included in any studies of this problem. 

We thank A. Cumming, K. Freeman, G. Novak, and P. Sackett for useful discussions. 
This work was supported by the Victorian Partnership for Advanced computing Expertise 
Program Grant Scheme, and by NASA through NAGS5-11779 under its Origins program, 
JPL 1228184 under its SIM program, and NSF through AST-9987417. 
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Fig. 1. — Eccentricity excitation by a companion planet, a, eccentricity evolution in the absence 
of tidal damping. M* = O.2M , M t = 1M©, M = lMj, m = 0.02 AU, a = 0.7 AU, e = 0.2, 
and periastra initially aligned. Solid curves: top=no perturbing accelerations; bottom=relativistic 
potential only. Dashed curves: top=spin and tidal bulge only; bottom=both relativity and spin 
and tidal bulges, b, eccentricity evolution with tidal damping: circulatory and libratory behaviour 
in the e^ — (vj\ — w ) plane. System parameters the same as in (a) for the case with no perturbing 
accelerations, but with initial e, = 0.01, periastra antialigned, and an unrealistic Qi = 0.2 used to 
illustrate the behaviour. In general, a two-planet system will evolve to a constant e^ with periastra 
aligned as long as r e is shorter than the lifetime of the system, c, dependence of e^ on Mi/M 
with 7 = 0. Solid curves: Newtonian point mass equations; dashed curves: Eqn. (4). Top set: 
a>o/ a i = 10; bottom set: a /ai = 20. d, dependence of e^ on Mi/M with relativistic effects 



included, a^/ai 
and (0.5,1). 



10. Bottom to top: (oi/AU, m;/Mj)=(0.05, 0.0033), (0.5, 0.0033), (0.05, 1) 



-11- 



00 

d 



< CD 
O 



o 



o 



d 



0.2 0.4 0. 



e 



oo 
d 



< CD 

d 

o 
D 



O 



d 




o 

II 
CD 
.1 



> 



Z5 

00 



o 



lo 



0.2 0.4 



0.6 



e 



Fig. 2. — Orbital decay timescales with and without relativity as functions of the outer 
eccentricity, e , and semimajor axis, a . M* = 1M , Mi = 1M®, M = lMj, a, = 0.03 AU. 
Top panel: no relativity. Bottom panel: with relativity. 
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Table 1: Structural parameters for various bodies 

Body fci Q Q' reference 



Venus 


0.295 


40 


203.4 


(1) 


Earth 


0.3 


21.5 


107.5 


(2) 


Mars 


0.14 


86 


921 


(3) 



Jupiter 0.34 2.2 x 10 5 10 6 (4) 

Sun 0.02 1.3 x 10 4 10 6 (5) 

M dwarf 0.30 2.0 x 10 5 10 6 (5) 

Love numbers (&/,), specific dissipation functions (Q), and modified Q- values (Q'), with 
Q' = 3Q/2kL (Goldreich & Soter, 1966, Murray & Dermott, 1999). Love numbers for 
gaseous bodies assume polytropic structure with polytropic indices 1, 3 and 1.5 for Jupiter, 
Sun and M dwarf respectively. References. (1) Yoder, 1995a, (2) Dickey et al, 1994, (3) 
Yoder, 1995b, (4) Yoder & Peale, 1981, (5) Mathieu, 1994. 
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Tablc 2: Orbital decay timescales 



System 


M^/Mq 


Mi/Mj 


p . 

J spin 


a t (AU) 


■* orb 


Pr 


T e 


Ta(0) 


Ta(0.01) 


Ta(0.1) 


Sun/Jupiter 


1.00 


1.00 


28 


0.02 


1.04 d 


0.12 


0.03 


0.16 


0.16 


0.06 










0.03 


1.91 d 


0.08 


0.45 


2.34 


2.30 


0.88 










0.04 


2.93 d 


0.03 


2.90 


15.8 


15.6 


5.83 


Sun/Earth 


1.00 


0.003 


28 


0.01 


8.8 hr 


0.23 


2.2-5 


0.53 


0.009 


8.1-5 










0.015 


16.2 hr 


0.15 


3.0-4 


7.41 


0.12 


0.001 










0.02 


25.0 hr 


0.12 


0.002 


48.7 


0.78 


0.007 










0.03 


1.9 d 


0.08 


0.027 


oc 


10.9 


0.10 










0.05 


4.1 d 


0.05 


0.75 


oc 


oc 


2.82 


M dwarf/Earth 


0.2 


0.003 


40 


0.007 


10.3 hr 


0.10 


2.4-5 


0.97 


0.01 


8.9-5 










0.01 


17.6 hr 


0.07 


2.4-4 


9.88 


0.097 


9.0-4 










0.012 


23.1 hr 


0.06 


7.9-4 


34.4 


0.32 


0.003 










0.02 


2.32 d 


0.03 


0.02 


oc 


8.82 


0.082 










0.04 


6.56 d 


0.01 


1.97 


oc 


oc 


7.39 


OGLE-TR-56 


1.04 


0.9 


_ 


0.023 


1.212 d 


0.11 


0.02 


0.40 


0.37 


0.05 


HD83443 


0.79 


0.41 


36.5 


0.038 


2.986 d 


0.05 


2.1 


73.2 






HD46375 


1.00 


0.25 


- 


0.041 


3.024 d 


0.06 


5.1 


26.7 






HD209458 


1.05 


0.66 


- 


0.05 


3.525 d 


0.05 


0.84 


59.8 


51.9 


3.43 



Orbital decay and eccentricity damping timescales (Gyr), and transit probabilities (Pr), for 
various real and hypothetical single planet systems. M^/Mj is the minimum mass of the 
planet in units of the mass of Jupiter, the actual mass being a factor (sini) -1 higher when 
% is the inclination of the orbit to our line of sight, 'd' stands for days, and the argument 
of r a is the eccentricity. For all stars we took Q'^ = 10 6 , for the known systems and the 
hypothetical Jupiter systems Q\ = 10 6 , while for the hypothetical Earth systems Q\ = 21.5. 
The symbol oo corresponds to r a > 10 11 yr. Entries in boldface for the last M dwarf-Earth 
system correspond to the approximate habitable zone for such a star. 



